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Summary 


The Jet Engine IR Analysis Procedure (J-EIRP) is a set of computer programs 
developed to evaluate jet engine cavity hot parts source radiation. These programs 
were intended for use on a personal computer, but include the ability to solve large 
problems on a mainframe or super-computer. The goal is to introduce the 
terminology and solution process used in J-EIRP, and provide insight into the radiation 
heat transfer principles used in this procedure. A sample jet engine cavity analysis 
demonstrates the procedure and capabilities within J-EIRP, and is compared to a 
simplified method for approximating cavity radiation. 


Introduction 


Thermal radiation from partially enclosed jet engine cavity surfaces can lead to 
complex analysis methods due to nonuniformly distributed surface radiant energy 
resulting from reflecting radiation and non-uniform surface cooling. Depending on 
engine internal surface geometry and thermophysical properties, the energy leaving a 
cavity can be composed of direct emission and reflected energy. Following the 
propagation of energy within the cavity becomes an enormous problem as the 
complexity of the geometry increases. Applying computer models and numerical 
approaches to evaluate the emission and propagation of multi-reflecting radiation can 
be an involved process. Currently, there are a variety of computer programs and 
simplifying analytical techniques used to analyze radiation as it propagates from a jet 
engine cavity. Each computer program has advantages and disadvantages based on 
the fundamentals of physics, modeling requirements and limitations, computer 
requirements, solution convergence requirements and computer limits. 


This paper describes an analysis procedure used in conjunction with a personal 
computer to evaluate the thermal radiation from complex jet engine cavities. The Jet 
Engine IR Analysis Procedure (J-EIRP) was developed to give the user flexibility in 
analyzing both simple and complex jet engine cavity geometries. A major objective of 
this effort was to develop a highly accurate thermal analysis procedure that can 
analyze a variety of radiation heat transfer problems with limited computer resources. 
J-EIRP was designed to manage components of radiation heat transfer problems on a 
personal computer, but also include the ability to solve more computer intense 


problems on a mainframe or super-computer. J-EIRF also provides the user a tool for 
developing thermal design experience and engineering judgment through analysis 
experimentation, while using minimal computer resources. 


Symbols 


A surface area 

B radiation interchange factor 

C V C 2 Planck's spectral energy distribution constants 
C Q speed of light in vacuum 

e emissive power 

E emissive energy per unit time 

F view factor 

h Planck's constant 

i radiant intensity per unit area 

I radiant intensity 

source radiation per unit solid angle ( PIREP notation) 
k Boltzmann constant 

n nth surface 

r radius distance 

T absolute temperature 

o Stefan- Boltzmann constant 

e emissivity 

6 polar angle 

4> circumferential angle 

T emissive power fraction in spectral region 
6 identity matrix 

k wavelength (vacuum) in micrometer (n/n) 

Superscript 

/ directional quantity 


Subscripts 

A of surface A 

b biackbody 

I, J ith and Jth surfaces 

k spectrally dependent wavelength 

« infinity 
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Method of Analysis 


The computer programs that comprise J-EIRP are displayed in flow chart form 
in figure 1 . The following brief descriptions of the computer codes exposes the 
required steps used in the J-EIRP solution process. Appendix A contains background 
information on radiation heat transfer terminology and concepts applied in the 
following programs. Appendix B contains the general solution method used in J-EIRP. 

Model 


The term "model" represents a data set that contains required information 
that defines the physical attributes of the jet engine cavity for analysis. A 
model supplies parameters that represent the physical jet engine cavity 
geometry, surface thermophysical properties, and solution requirements. 
The cavity geometry is subdivided into a number of finite surfaces called 
nodes. The actual number of nodes and their geometry depends on the 
desired model representation, accuracy of results, structural design 
considerations, computer capabilities, and computer computational time 
requirements. Each node represents an average surface temperature 
defined over the entire surface. Using the ray tracing approach, 
restrictions on uniform node energy distributions (radiosity and 
irradiation) may be eliminated (ref.1). 


NEVADA 

The Net Energy Verification and Determination Analyzer program or "NEVADA" 
(ref. 2), was selected to simulate the radiation propagating from modeled 
surfaces. NEVADA is a software package consisting of several programs in 
which a Monte-Carlo mathematical technique is applied to radiation 
propagation. The NEVADA program was attractive to use because of its ability 
to model the complex laws of physics associated with radiation propagating 
from various geometry defining surface types and thermophysical properties. In 
the NEVADA code, a statistical numerical method using the Monte-Carlo 
technique is applied to a ray tracing procedure to model radiation exchange. 
The ray tracing procedure mathematically traces emitted rays (simulating 
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emitted radiation) as they propagate. Each ray leaving a surface is considered 
a bundle of photons which carries equal, discrete, amounts of energy. The 
path of the bundles (rays) may interact with various surfaces, which may 
reduce the energy level of the bundle. The interacting surfaces may have 
different thermophysical properties and geometric configurations that may affect 
the propagation of each bundle differently. By accounting for all the emitted 
bundles as they propagate, percentages of incident and absorbed energy at 
different locations can be computed. The percentages of absorbed energies 
(radiation interchange factors) are then applied to the energy balance 
equations. 


NEVADA is an industry standard code used in predicting radiation interchange 
factors. Various documented uses and results are published. The NEVADA 
numerical ray tracing technique as applied to cavity radiation propagation for 
isothermal and nonisothermal gray-diffuse cylindrical geometries with uniform 
surface emissivity values was verified in reference 1. 


TRACK-E 

The Track Energy Procedure or "TRACK-E“ is a FORTRAN program used to 
evaluate the energy propagation between defined surface nodes. The TRACK- 
E program calculates and sums the surface energy emitted and absorbed 
through the radiation conductor values (term used to represent the energy flow 
paths between separated surfaces). The radiation conductors incorporate the 
NEVADA radiation interchange factors. The TRACK-E program can also modify 
the radiation conductor values to include wavelength band evaluation of the IR 
spectrum. If the NEVADA model includes spherical or hemispherical surfaces 
that enclose the cavity opening for energy propagation studies, TRACK-E can 
then be used to evaluate the total and component energy distribution over 
these surfaces. Appendix B defines the general radiation exchange equation 
solved in the TRACK-E program with the radiation interchange factors from the 
NEVADA program. The user can also evaluate radiant energy calculations for 
individual surfaces by entering radiation interchange factors, areas, emissivities, 
and temperatures. The program also tabulates the emitted and absorbed 
energy from each surface. 


Sample Calculation 


The use of J-EIRP to evaluate the radiation emission from a conceptual jet 
engine cavity will now be described. This example problem will serve to illustrate the 
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programs and computer requirements to achieve radiation wavelength band solutions. 
This sample cavity analysis will model diffusely reflecting radiation propagating from a 
jet engine cavity to its surrounding environment using a personal computer (80386 
microprocessor, 33 MHz speed) for all calculations. Figure 2 displays a computer- 
generated representation of jet engine cavity surface geometry. This jet engine 
surface geometry is defined by eight surface nodes. These eight nodes accurately 
represent the jet engine cavity because of the cavity symmetry and uniform surface 
temperatures over the individual node regions. An 18-node hemisphere, spherically 
sectioned, and placed behind the cavity exit plane is used to absorb and evaluate the 
incident radiant energy from the cavity surfaces. The computer-generated 
representation of the hemisphere is shown in figure 3. This sample cavity analysis will 
display several J-EIRP features. Various other options exist within J-EIRP, including a 
variety of ray tracing options within NEVADA, which may be beneficial depending on 
the scope and requirements of the analysis. 


The NEVADA ray tracing analysis computes the view factors and radiation 
interchange factors between the cavity and hemispherical surfaces. These results 
relate the percent of energy emitted and reflected between the cavity surfaces, and the 
cavity and hemispherical surfaces. Figure 4 is a polar plot displaying the radiation 
interchange factor results for two jet engine cavity components to the hemispherical 
surfaces. The solid line represents radiation from the jet engine turbine geometry 
while the dotted line represents radiation from the jet engine divergent nozzle. These 
results show the relationship of component geometry and surface properties on the 
engine radiation emission before introducing surface temperatures. These results may 
be useful for evaluating the influence of geometry and surface properties on the 
transfer of energy to the hemispherical sections. In this example, the 90 degree 
location represents the axial direction or angle normal to the jet engine cavity exit 
plane. Due to the statistical nature of the ray tracing analysis, computer computational 
time depends on the desired accuracy. These results represent solutions with small 
error bands (high accuracy), which required roughly two hours of computer 
computational time. Useful but less accurate results require much less computer 
computational time. 


The TRACK-E program now allows one to assign temperatures to the engine 
cavity surfaces and evaluate the radiant energy emission from the total jet engine 
cavity, groups of cavity components, or individual cavity components. The user has 
the option to evaluate a wavelength band of the IR spectrum or the total IR energy 
transfer. In this example the TRACK-E program will be used to evaluate several 
wavelength bands and the total IR energy transferred to the hemisphere. The TRACK- 
E program reads the NEVADA output data set that contains surface areas and the 
radiation interchange factors. A second file providing surface temperatures is required 
for evaluating the total or wavelength band energy transferred between the surfaces. 
The user has the option of reviewing J-EIRP results in terms of energy through each 
radiation conductor, net energy results for individual surfaces, total hemispherical 
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energy results, total cavity energy results, and tables of hemispherical results used for 
graphing. Polar plots can also display the radiant energy emission results in the form 
of watts per unit solid angle (steradian). Figure 5 displays the total radiant energy 
emitted by the cavity and figure 6 displays two wavelength band results of interest. 
Figure 7 compares the divergent nozzle component energy emission to that of the 
total jet engine cavity (radiation interchange factors for this case were shown in fig. 4). 
The importance of component analysis is seen in figure 7. The divergent nozzle 
contributes nearly all the off-axis radiation (0.0 to 45.0 and 135.0 to 180.0 degrees). 
This sample TRACK-E cavity analysis, which evaluated 160 NEVADA radiation 
interchange factors, required less then five minutes of computer computational time. 


A simplified jet engine hot parts radiation analysis method called PIREP (ref.3), 
described in appendix C, can now be compared to the J-EIRP results for the direction 
normal to the exit plane. The same jet engine configuration was analyzed using the 
PIREP method for the 3.0 to 5.0 micrometer wavelength band. This PIREP analysis 
resulted in 204.6 watts per steradian. This value corresponds to the J-EIRP results of 
291.7 watts per steradian, as shown in figure 6. The PIREP procedure does not have 
the ability to calculate off-normal radiation emission or individual engine component 
contributions to the total radiation emission. 


Concluding Remarks 


A thermal radiation analysis method called Jet Engine IR Analysis Procedure (J- 
EIRP) was developed to evaluate jet engine cavity hot parts source radiation. The 
objectives of the J-EIRP development were to achieve the greatest accuracy in model 
representation and solution, while minimizing computer resources and computational 
time. The computer programs that comprise J-EIRP were selected on the basis of 
their performance, accuracy and flexibility to solve both simple and complex problems, 
while retaining the ability to add or substitute programs based on user requirements. 
These programs were intended for use on a personal computer, but include the ability 
to solve large problems on a mainframe or super-computer. J-EIRP was also 
designed so programs could be added or substituted based on user requirements. 
These options produce a working environment for solving complex problems that 
require limited or minimal computer resources. J-EIRP also provides the user a tool 
for developing thermal design experience and engineering judgment through analysis 
experimentation. 


6 



Appendix A - Basic Radiation Heat Transfer Concepts 


The following discussion briefly introduces some basic radiation heat transfer 
concepts. These radiation heat transfer concepts are used to define the solution 
technique and requirements for use in J-EIRP. The basic definitions, equations, units, 
and notation used to define the radiation characteristics were selected from reference 
4. 


The simplest mode of radiation heat transfer occurs when energy is emitted 
from a single flat surface in a vacuum where no reflected energy is involved. If the 
surface is a perfect absorber and emitter of radiant energy, it is called a blackbody. 
This blackbody also emits the maximum amount of radiant energy to its surrounding 
hemisphere. The intensity of radiation emitted by a blackbody in any direction for a 
single wavelength is defined by Planck’s Law as: 




gg, 

X 5 (e^ r -1) 


BTU 

hr \im ft 2 sr 


(A1) 


where C, and C 2 represent Planck’s spectral energy distribution constants, as shown: 


C , « hC * 0 



hCo 

k 


C. = .18878X/0 8 gTU 

hr ft 2 


C 2 = 25897.84 ji m°R 


The prime notation defines a directional quantity of radiation per unit solid angle for a 
single direction. The subscripted terms specify the quantity as spectral (X - one 
wavelength), and the surface as a blackbody (b). The dependent variables are listed 
in the parentheses. The emissive power radiated from a blackbody over an entire 
hemisphere in a vacuum at a particular wavelength is defined as: 


eJKT) - vi’JX.T) = 


2n C, 

X 5( d cyir_i) 


BTU 

hr \xm ft 2 


(A2) 


By integrating equation (A2) over the entire wavelength band, the total blackbody 
hemispherical emissive power at a specific surface temperature is then: 
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(A3) 


9 b { T) - J &xb^ A) ~ f 11 “ f 


2nG 


o JL B (a <v * r -1) 


-<*W = 0 7^ 


arc/ 

hr ft 2 


where o represents the Stefan-Boltzmann constant: 


a = 0.17123x70-® 


BTU 

hr ft 2 °/? 4 . 


Applying the emitting surface area to equation (A3) gives the total blackbody 
hemispherical emitted energy per unit time as: 


- f 0 Ai&xbM^ ~ 


B7U 

hr 


(A4) 


For a real (non-blackbody) surface, thermophysical properties control the 
radiant energy absorbed, emitted, and reflected from the surface. These properties 
can be a function of incident and reflected angles, wavelength, and temperature. 
Frequently, several simplifying assumptions can be made regarding the surface 
thermophysical properties. Averaged surface properties may be derived if one can 
average over all directions and wavelength bands. Applying these assumptions, a 
surface that absorbs and emits a fixed fraction of radiation from any direction and at 
any wavelength is now defined as a diffuse gray surface. The directional and spectral 
absorptivity and emissivity then become: 

«'x(B,4>XTJ = a A (TJ 


From Kirchhoff’s law the directional, spectral and hemispherical total values of 
absorptivity and emissivity are equal, thus: 

e(TJ = 


(A5) 


The dependence of the surface temperature on the surface emissivity value also may 
be averaged over the radiation heat transfer band. 


Applying surface thermophysical properties results in the surface being defined 
as a non-blackbody. The total non-blackbody (diffuse-gray) hemispherical emissive 
power is: 
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(AS) 


^a oT a 


BW 
hr ft 2 


The total non-blackbody hemispherical emitted energy is: 

^A ,(7*) = f Q ^1 e A 0 Xb(^)^ = ^i e A°^A 


BTU 

hr 


(A7) 


For radiation to be exchanged between surfaces, a geometric relationship 
(geometric configuration factor or view factor) between the surfaces must be defined. 
The view factor for blackbody isothermal surfaces is the fraction of emitted radiant 
energy leaving a surface i that reaches another surface j. Since both surfaces are 
blackbodies or perfect absorbers of radiant energy, no reflected energy is introduced 
into the view factor value, F h ,. The total blackbody emitted energy from surface A, 
that is incident and absorbed by surface A 2 is then: 


a 1 - a£J/) - A Fj^ „ ^o(7^) 


BTU 

hr 


(A8) 


For simple geometric configurations various view factor references are available. For 
complex geometric configurations, various computer techniques have been developed 
where the surfaces are defined, and mathematical or ray tracing techniques are 
applied to solve the surface relationships. 


For non-blackbody surfaces the blackbody view factor is replaced by the 
radiation interchange factor, B hj , which represents real surface radiation exchange. 

The radiation interchange factor is the fraction of emitted energy by a real surface i 
that is absorbed by real surface j, including all reflections from other real surfaces 
including the emitting surface i. Computer programs are also available to calculate B H 
values. For the radiation exchange between two real surfaces, defined as A, and A 2 , 
the total non-blackbody emitted energy from surface A, that is absorbed by surface A 2 
is: 


^ 4 , - A^a) ~ ^ 1 ^ 4 , - A^Ay°^X 


BTU 

hr 


(A9) 


This is where the B hj term defines the relationship between the energy emitted from 
surface A 1 that is absorbed by surface A 2 directly and indirectly. This indirect radiant 
energy may be in the form of multi-reflecting energy between A 1 and A 2 , and energy 
from A 1 that is reflected from other surfaces to A^ The radiation path for a multi- 
surface cavity configuration may include all surfaces, which can result in a large 
number of total non-blackbody emitted energy terms: 
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Energy Transfer From All Possible Radiation Paths = £ Y G /°A®/-/(7?) 

M /=i 


The radiation exchange also can be evaluated within a particular wavelength 
region or fraction of the total emissivity power. This is done by simply using the total 
blackbody hemispherical emissive power, equation (A3), and modifying the limits of 
integration as shown: 


e b (XJ^k z 7) = Je xb d(X) = f 


2 nC, 


1, xV^-i) 


<*X) 


(A 10) 


The non-blackbody hemispherical emitted energy within a wavelength band becomes: 


^(*i^A 2 7) - e^A A J e xb d{X) = e^A A J 


2nC^ 


k s (e CJXT - 1) 


-4M 


(AH) 


Applying the radiation interchange factor results in the non-blackbody emitted energy 
from individual surfaces within a wavelength band: 




2nC< 


i, aV^-i) 




(A 12) 


The integration term can be solved using a numerical polynomial curve fit resulting in 
the total non-blackbody hemispherical emitted energy equation as: 




STL/ 


/?r 


(A13) 


where the r aAT is represented (ref. 4) as: 

15 


o-xr 


n n-1 


' & -nx 



n 

k n n z n 3 j 



(A14) 


w/we x = — 
AT 


Tabulated values of these polynomial approximations are also available as a function 
of wavelength-temperature products in reference 4. 
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Appendix B - The J-EIRP Solution Concept 


The programs that currently comprise J-EIRP solve the radiation heat transfer 
concepts described in appendix A. The same definitions, notations, and equations 
defined in appendix A apply, but a conversion factor modifying the units to watts has 
been included. The total non-blackbody hemispherical emitted energy in equation 
(A13) is the general form of J-EIRP solution concept. 


To absorb, identify, and evaluate the radiation from an engine cavity, a 
sectioned hemisphere placed behind the cavity exit plane can be used. This multi- 
sectioned hemisphere reveals the uniform intensity over the finite hemispherical 
surfaces areas from the jet engine cavity. By selecting a distance where the emitted 
and reflected cavity radiation results in nearly normal intersections with the 
hemispherical surfaces, the radiant energy can be calculated as a function of solid 
angle. The equation describing the radiant energy emitted from a defined cavity 
surface to a defined hemispherical surface in terms of a solid angle is as shown: 




2 ^ 


]o72, (.2931) 


f'VJ 


Watts 

sr 


(B1) 


In this equation r represents the distance from the source cavity surface area A, to the 
target hemispherical surface area A 2 for a solid angle calculation. By calculating the 
radiant energy emitted and reflected from each jet engine cavity surface to each 
hemispherical surface and summing the hemispherical results, the jet engine cavity 
radiation emission is calculated. Equation (B1) summarizes all the basic concepts 
within J-EIRP. The radiation interchange factors (B Hj values) are calculated using the 
NEVADA program which analyzes a model representing the jet engine cavity. The 
summed radiant energy from all the jet engine cavity surfaces to each hemispherical 
surface is calculated using the TRACK-E program. TRACK-E also evaluates the 
component radiant energy distributions over the hemispherical surfaces, and the 
radiant energy emitted and absorbed within the jet engine cavity. 
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Appendix C - The PIREP Solution Concept 


The Preliminary Infrared Radiation Emissions Program (ref. 3) includes a 
simplified method of calculating the jet engine cavity hot parts emission. The general 
PIREP method applies Planck’s Law directly to a solution where the jet engine cavity 
geometry is simplified to represent a flat plate emitting energy at the size of the 
exhaust nozzle throat. Therefore, this method evaluates only the directly emitted 
energy from a flat plate and eliminates cavity wall reflections. To better understand the 
general PIREP equation and results, the PIREP jet engine hot parts analysis method 
will be derived in the following discussion. 


The intensity of radiation emitted in any direction for a single wavelength is 
defined by Planck’s Law in equation (A1) as: 




x 5 (e^ r -1) 


BTU 

hr \im ft 2 sr 


Applying Planck’s spectral energy distribution constants as given in appendix A and 
converting units, equation (A1) becomes: 


I'uJLW 


2(0.18878x70®) 

.2931 ) 

j 76849.1917 

Watts 

A 5( e 25897.84/X7_ 1) 

144 J 

1 ‘ ^5^25897.84 

\im in 2 sr 


(Cl) 


This is the general form of Planck’s Law used in PIREP. Planck’s Law, which is for a 
single wavelength calculation, is normally integrated over the IR spectrum or 
wavelength band. To simplify the wavelength band integration, the integration is 
replaced with a wavelength delta. This assumes the wavelength band is small and 
that an average wavelength can be substituted into Planck’s Law as follows: 
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V i\ T i T\ - f 76849.1917 ^ _ 76849.1917 AX 

“ 2 ' j ^5^25897.84/i7_^ ~ j5(e 25897M A 7 '- < |) 

where: 

AX = wavelength difference 
X = average wavelength 


<C2) 


This represents the intensity of emitted radiation per unit solid angle from a flat plate in 
all directions, where temperature is defined by the turbine exit total temperature. By 
introducing the exhaust nozzle throat area into equation (C2), PIREP evaluates the 
intensity from a jet engine cavity. This is shown in the following equation along with 
the notation changed to reflect PIREP’s form. 


// i\ t \ t\ - J *° - 76849.1917 AX 
bAr ' 1 ~ 2 ' A s ~ x B (e 25 * 97M *T-V 


Watts 
in 2 sr 


(C3) 


The PIREP method was designed as a fast approximation tool for calculating the 
emission from jet engine cavity hot parts. Alternative solution methods are required if 
component effects are to be evaluated, or increased analytical accuracy is required. 
Equation (C3) differs slightly from the actual PIREP equation due to the use of updated 
values for Planck’s spectral energy distribution constants (C, and C 2 ). 
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SUPPLY SURFACE 
TEMPERATURES 


Figure 1. - Jet Engine IR Analysis Procedure (J-EIRP). 




Figure 2. - Surface geometry of a sample cavity model. 
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Figure 3. - Surface geometry of sample hemisphere model section 
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Figure 4. - Sample cavity component radiation interchange factors. 
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Figure 5. - Sample cavity total radiation emission. 
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Figure 6. - Sample cavity wavelength band radiation emission. 
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Figure 7. - Sample cavity divergent nozzle effect on radiation emission. 

(3.0 ^m< wavelength <5.0 urn) 
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